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ABSTRACT 

X-ray spectra in the range 1.5 — 8.5 keV have been analyzed for 526 large 
flares detected with the Solar Assembly for X-rays (SAX) on the Mercury MES¬ 
SENGER spacecraft between 2007 and 2013. For each flare, the temperature and 
emission measure of the emitting plasma were determined from the spectrum of 
the continuum. In addition, with the SAX energy resolution of 0.6 keV (FWHM) 
at 6 keV, the intensities of the clearly resolved Fe-line complex at 6.7 keV and 
the Ca-line complex at 3.9 keV were determined, along with those of unresolved 
line complexes from S, Si, and Ar at lower energies. Comparisons of these line 
intensities with theoretical spectra allow the abundances of these elements rela¬ 
tive to hydrogen to be derived, with uncertainties due to instrument calibration 
and the unknown temperature distribution of the emitting plasma. While signif¬ 
icant deviations are found for the abundances of Fe and Ca from flare to flare, 
the abundances averaged over all flares are found to be enhanced over photo- 
spheric values by factors of 1.66 ± 0.34 (Fe), 3.89 ± 0.76 (Ca), 1.23 ± 0.45 
(S), 1.64 ± 0.66 (Si), and 2.48 ± 0.90 (Ar). These factors differ from previous 
reported values for Fe and Si at least. They suggest a more complex relation of 
abundance enhancement with the first ionization potential (FIP) of the element 
than previously considered, with the possibility that fractionation occurs in flares 
for elements with a FIP of less than ~7 eV rather than ~ 10 eV. 

Subject headings: Sun: abundances — Sun: corona — Sun: flares — Sun: X-rays, 
gamma rays — line: identification 
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1. INTRODUCTION 


The Solar Assem bly for X-rays (SAX) is part of the X-ray Spectrometer (XRS: 


Schlemm et a 1. 


( 2007 1) on the MESS ENGER (MErcury Surface, Space ENvironment. 


(Solomon et al. 

2001; 

^ ? — r ^ 

Gold et al. 

2001 


Santo et ah 


2001 ). Its primary role is to measure the spectrum of solar soft X-rays that are the main 


exciters of fluorescence radiation from Mercury’s surface. The spectrum of this fluorescence 
radiation is measured, in turn, with three gas proportional counters of XRS to provide 


information on the chemical composition of the surface (INittler et al. 


20111). SAX consists 


of a cooled, 500 /ini thick Si-PIN solid-state detector covering the energy range from ~1 to 
~10 keV with an energy resolution of 598 eV at 5.9 keV. Its pinhole aperture has an area 
of only 0.03 mm 2 so that the detector count rates are limited to manageable levels during 
all but the most intense flares. The detector assembly is mounted on the MESSENGER 
sunshade and views the Sun directly through three Be foils (each 25/irn thick) to reduce the 
solar thermal input to the detector and attenuate the flux of the lowest energy X-rays. 


Since the launch of MESSENGER on 2004 August 3 and capture into orbit around 
Mercury in 2011, SAX has observed hundreds of solar flares and obtained soft X-ray 
spectra with a cadence generally of 5 minutes. The emission in the 1-10 keV energy range 
is dominated by free-free and free-bound thermal continuum radiation but also includes 
spectral line complexes that can be resolved with the SAX spectral resolution. These 
include complexes of lines due to highly ionized Fe (peaking at ~ 6.7 keV and ~ 8 keV) 
and Ca (peaking at ~ 3.9 keV). Estimates of the abundances of Fe and Ca can be made 
from spectral fits to the measured count-rate spectra with a precision depending primarily 
on the accuracy of the assumed temperature distribution. Forms for the distribution (the 
differential emission measure) which are a power-law or exponential in temperature were 
considered, but most of the analysis was done with simpler forms described by a single 
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temperature (isothermal) or two temperatures. At energies of < 3.5 keV, contributions to 
the emission are also made by lines of S, Si, Ar, Mg, and Fe ions but these are not resolved 
with SAX. Consequently, estimates of the S, Si, and Ar abundances can also be made but 
with greater uncertainty because of their smaller contributions to the total emission. Mg 
makes too small a contribution to the spectrum and at too low an energy for a meaningful 
abundance determination to be made from SAX spectra. This paper is primarily concerned 
with deriving estimates of the abundances of these elements for flares observed with SAX 
and comparison with those from other instruments and with theory. 


Element abundances in coronal plasmas such as flares are widely considered to be 
different from those in the photosphere by amounts that depend on the first ionization 
potential (FIP) of the element - the so-called FIP effect - but there is still debate about the 
magnitude of the difference for each element or the so-called FIP bias, equal to t he ratio 


of the measured element abundance to the abundance in the photosphere. 


Feldman et ah 


(11992 1 give values of element abundances in the corona that are as much as a factor of 
4 higher than photospheric abundances (a FIP bias of ~ 4) for elements with relatively 
low FIP (taken to be < 10 eV) but are approximately equal to photospheric abundances 


(FIP bias of 1) for high-FIP (> 10 eV) elements. The later work of 


F eldm an fe L am ing 


(j2000Q_ gives reduced abundances of some low-FIP elements, notably Mg, Si, and Ca. 


The 


Feldman et a! 


(119921) values are the basis for the “coro nal” abundance set in the 


CHIANTI database and so 


t ware p ackage used extensively here (IDere et al. 


1997; 


Landi et al. 


2013). Fludra & Schmelz ( 19991 ) propose coronal abundances (the “hybrid” abundance 


set in chianti) that are enhanced over the photospheric abundance estimates available 


at that time for low-FIP elements by 


for high-FIP elements. 


Schmelz et al. 


actor s of ~2 but reduced by factors of up to ~2 


I 2012) have given recommended abundances for the 


solar corona but this is an amalgam of results from coronal spectroscopy, the solar wind, 
and solar energetic particles so it may not be appropriate for solar flares. The most recent 
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photospheric abundance estimates (lAsplund et al.l 12009k ICaffau et al.l 1201 ll) have resulted 
in some minor revisions of elements considered here (Fe, Ca, S, and Si) compared with 
previous work (although there are substantial reductions in the C, N, and O abundances). 


The ph otosp heri c a 


from 


n undan ces used here to evaluate FIP bias values given in Table Q] are 


Asplund et ah 


(2009) 


More refined estimates of coronal abundances of S, Si, Ar, and K have been made 
from flare and non-flare o bserva tions of X-ray lines by the RESIK crystal spectrometer 
on Coronas-F (see 


Svlwes te r et al 


2014 and references the rein) and of Fe from RH ESSI 


observations of the line complexes at 6.7 keV and 8 keV .([Phil 


Compared with the photospheric abundances of 


Asplund et ah 


ips fe Dennis 


I 


2012) 


(120091 ). they indicate FIP 


biases which are listed in Table [H Th ere is some disagreement betwe e n the se estimates 


and those of 


Fludra & Schmelz (119991 ) (for S) and 


Feldman & Laming (2000) (for Fe). The 


RESIK and RHESSI sp ectra t hus sugge st t hat the dependence on FIP is not as simple as 


that given originally by 


Feldman et al. 


(119921 ). This is also indicated by some more recent 


observational results which we will discuss later. 


Theoretical models to explain the FIP effect generally use the fact that low-FIP 


elements exist in a partially or wholly ionized form in the photosp 
but high-FIP elements are neutral. Earlier models, summarized by 


Here or c. 


Henoux 


iromosphere 


199a), include 


mechanisms like diffusion, with magnetic fields playing an active or passive role. In the 


model developed by 


Laming (2012) (and references therein), FlP-dependent abundance 


enhancements occur as the result of ponderomotive forces associated with Alfven waves 
excited in the solar corona, and can be calculated for different values of coronal flux loop 
lengths and wave flux. Thus, comparison with observations is possible in principle. 


In this paper, we report on an extensive analysis of MESSENGER SAX spectra taken 
during 526 well-observed flares detected between 2007 and 2013, when the spacecraft was 
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first in its cruise phase to Mercury and then in orbit about Mercury. Abundance estimates 
for the elements Fe, Ca, S, Si, and Ar are given and their distributions examined. The 
observations and analysis procedures are described in Section 2 for both a single flare 
and for all flares in the data set with sufficient intensity to detect emission in the Fe-line 


complex at 6.7 keV. The resulting abundance values are given in Section 3, anc 


in Section 4 with previous results and with the predictions of the 


Laming; (2012) theory. 


compared 


2. OBSERVATIONS AND DATA ANALYSIS 
2.1. Flare Selection and Time Coverage 

SAX data were normally collected continuously over 5-minute intervals although on 
some occasions the instrument was commanded into a high-cadence mode with 20-s or 40-s 
data-gathering intervals close to periapsis of the orbit about Mercury and during some of 
the more intense flares. Only those time intervals in which the count rate was > 5 000 s” 1 in 
the 6.3-7.0 keV energy range were selected for analysis to ensure that the Fe line complex at 
6.7 keV was well observed. This resulted in the selection of 656 flares or significant increases 
in X-ray flux for analysis. Of these, events lasting for only a single 5-minute data-gathering 
interval were eliminated from this data set leaving a total of 526 flares with at least two 
contiguous data intervals that were used to determine element abundances. The results 
from these flares were used in computing the average abundances listed in Table [2] A full 
catalog of flares analyzed with best-fit spectral parameters and more extensive plots are 
available at the following location: 

http: //hesperia.gsfc.nasa.gov/messenger/results/plots_and_tables/. 

The analyzed flares covered the period from 2007 May 28 to 2013 August 19. They 
are not uniformly distributed in time for three main reasons: the varying solar activity 
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through the deep minimum in 2009, intermittent detector operations necessitated by the 
various spacecraft maneuvers and periods of high detector temperatures, and particle 
contamination when the spacecraft was presumably well connected magnetically to the 
flaring active region. As a result, there are 27 flares analyzed up until 2007 December 31 


and the remainder after 20 
been given by Ni ttler et al 


0 February 6. An analysis of eleven flares in 2011 has previously 


(201l|). Four of these satisfy our selection criteria and so are 


included in our lists of analyzed events. The others were either too small to be included or 
occurred during a period of high particle contamination. 


2.2. Analysis of Flare on 2007 June 1 

A detailed analysis of an M2.8 flare observed by SAX on 2007 June 1 is first given 
as an illustration of the methods we used. Figure [1] shows X-ray light curves as recorded 
by SAX and GOES. Since the event was from NOAA active region 10960 at S08E78 and 
Mercury was at a Stonyhurst longitude of 49° East, both instruments would have seen the 
event equally well. The GOES light curves show an initial peak at 14:22 UT followed by 
the main event peaking at 14:59 UT, then a gradual decay to the pre-flare level at about 
16:15 UT. As would be expected from their sensitivity to similar energy ranges, the SAX 
light curves (corrected for the light travel time from MESSENGER to 1 AU) closely match 
the GOES light curves, although with a time resolution of 300 s compared with 3 s for 
GOES (2 s after 2009 December 2 with GOES 14 and GOES 15). 

We describe the fitting procedure for the spectrum taken during the 5-minute interval 
beginning at 14:56:29.9 UT, at flare maximum. The observed spectrum, shown in the 
two panels of Figure [21 is dominated by continuum (the sum of free-free and free-bound 
emission), with line complexes due to Ca ions (~ 3.9 keV) and Fe ions (~ 6.7 keV and 
~ 8 keV) also evident. The temperature governs the slope of the spectrum when plotted 
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Fig. 1.— SAX and GOES light curves for the M2.8 flare on 2007 June 1. Top panel: SAX 
light curves plotted with 300 s time bins in four energy bands: 1.3 — 2 keV (black), 2 — 3.2 keV 
(pink), 3.2 —5.3 keV (green), 5.3 —8.5 keV (cyan). Bottom panel: GOES X-ray light curves 
at 1 — 8 A (black) and 0.5 — 4 A (pink) with 3 s time bins. 
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against energy, with higher temperatures giving rise to spectra with flatter slopes, while the 
emission measure determines the absolute value of the flux at any energy. 


The measured count rate spectrum was fitte d wit 


from the CHIANTI package (v. 7.0) (jLandi et ah 


1 theoretical photon spectra calculated 


2013). These spectra are defined by values 


of the electron temperature, T, and the volume emission measure, EM = N^V (where 
N e = electron density, V = emitting volume), for a set of abunda nces w hich in our analysis 


were chosen to be relative to t 
fractions from 


Mazzotta et ah 


re coronal abundances of 


Feldman et al 


(1992) with ion 


( 19981 ). In order to speed up the calculation of the theoretical 


thermal spectra, base isothermal spectra were obtained from an abstraction of the CHIANTI 
routines. Tables were generated for each chianti spectral line for a logarithmically spaced 
temperature grid with the grid spacing dependent on the significance of that line. When a 
particular spectrum was generated from that database, the lines were found in the energy 
range of interest and the emissivity computed by interpolating from the table. The total 
line spectrum was then computed by summing over all line contributions with each line 
weighted by its assumed clement abundance. The continuum was similarly computed and 
made available either as a separate spectrum or combined with the total line spectrum. In 
this way, EM, T, and the abundances of Fe, Ca, S, Si, Ar, and other elements could all be 
used as free parameters and adjusted as necessary in the fitting procedure. 

As the temperature structure of the flare emission is unknown, functional forms must 
be assumed. The simplest case is an isothermal function defined only by T and EM. 

This appears to apply most nearly to X-ray spectra over limited energy ranges taken 
near maximum and in the decay of a flare, as has been shown in analysis of GOES flares 


by 


Garcia! (119941) and with tem perature-dependent s pectral line ratios fr om the Yohko h 


Bragg Crystal Spectrometer by 


Phillips et al. 


( 20051 ). In the analysis of 


Niftier et ah 


(20111). SAX spectra were fitted assuming an isothermal emitting region for the 1 — 8 keV 
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spectral range on the grounds that two-temperature models gave worse fits. However, 
isothermal fits to the entire spectral range of SAX are not likely to be valid in principle 
for most spectra. This is because previous flare observations with high-resolution crystal 
spectrometers have shown that temperatures determined from the Ca lines at 3.9 keV (due 
to Ca XIX lines with relatively weaker satellites) are found to be lower by typically ~ 5 MK 
than those determined from the Fe-linc complex at 6.7 keV (duetto Fe XXV lines with 


temperature-dependent dielectronic satellites) (e.g. 


Feldman et al. 


1980). An isothermal 


assumption is more likely to be valid over limited spectral regions, particularly on either side 
of the 6.7 keV Fe-line complex, where the neighboring continuum has a slope determined by 
a temperature which is likely to approximate that derived from satellite-to-resonance line 
ratios in the Fe XXV lines and nearby satellite lines. 

For the SAX spectrum in Figure [21 four forms of the temperature dependence of the 
differential emission measure, defined as DEM = N^dV/dT, were assumed: (a) an isothermal 
model (“1-T”) described by single values of T and EM; ( b ) a two-temperature model (“2-T”), 
with two values of T and EM; (c) a multi-thermal model with power-law dependence 
of the differential emission measure on T (“multitherm_pow”), constant xT~ a (where 
a is positive); (d) a multi-thermal model with the differential emission measure having 
an exponentially decreasing dependence on T (“multitherm_exp”), constant xexp (—T/r) 
(where r is positive). In our fitting procedure, theoretical spectra calculated from chianti 
were folded through each of these functional forms with free parameters (e.g. T, EM, a, and 
r) and the abundances of elements significantly affecting the spectrum in this energy range, 
viz., Fe, Ca, S, Si, and Ar. For the 2-T model the abundances of the higher temperature 
component were allowed to vary but the abunda nces o f the lower temperature component 


were held fixed at the chianti “coronal” values ( Feldman et al 


1992) listed in Table [D 


The fitting routines are much the same as those used by 


Phillips & Dennis (2012) to ob- 
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tain the Fe abundance from RHESSI spectra. Data files containing the count-rate spectrum 
and detector response matrix were read by a program called Object Spectral Analysis Execu¬ 
tive (OSPEX) (http://hesperia.gsfc.nasa.gov/rhessi3/software/spectroscopy/spectral-analysis-software/). 
Theoretical spectra calculated from each of the four model spectral functions (a-d) were 
then folded through the detector response function with first-guess values of the free 
parameters (e.g. for model (a) T, EM, and abundances of Fe, Ca, S, Si, and Ar) to calculate 
predicted count rates. Fits to the background-subtracted SAX count-rate spectrum (with 
the background determined before or after the flare) were accomplished with OSPEX by 
continual adjustment of the fitting parameters to improve the goodness of fit, determined 
by the reduced chi-squared statistic - 


2 _ 1 \ 

Xred ., 


V 




[(C, 


i.obs 


rv \2 / 2 
i,exp) / 


i,expi 


(i) 


where u is the number of degrees of freedom (number of data points minus the number of 
free parameters), Ci )0 b s and are the measured and predicted count rates respectively in 
the ?’th energy bin, and a ijSxp is the standard deviation of the expected count rate assuming 
Poisson statistics. 


Figure [2] shows the 1-T and 2-T fits to the measured SAX count rate flux spectrum with 
normalized residuals below. (Equivalent plots for the multitherm_pow and multitherm_exp 
fits are not shown but have similar appearance and all have acceptable values of xf ed .) 
Count rate flux is the measured count rate converted to units of counts s -1 cnr 2 keV -1 by 
dividing by the effective area projected to 1 AU and by the energy bin width. This enables 
the varying distance of the MESSENGER spacecraft from the Sun to be accounted for as it 
passed from Earth to Mercury and during its elliptical orbit about Mercury. For the 1-T 
fit, the energy range (indicated in the figure) was confined to 2.3 — 8.5 keV because a good 
fit could not be achieved at lower energies. For the 2-T function, a satisfactory fit was 


1 Temperature 


2 Temperature 





Energy (keV) 


6 7 8 9 10 





Energy (keV) 


6 7 8 9 10 


Fig. 2.— SAX background-subtracted count flux spectra (black histograms) and two model 
fits for the 14:56:30-15:01:29 UT period at the maximum of the 2007 June 1 flare with 
normalized residuals (the differences between the measured and best-fit fluxes as a fraction 
of the statistical la uncertainties) shown below each spectrum. Left: 1-T (isothermal) 
function (green curve) fitted over the 2.3 — 8.5 keV energy range with T — 22 MK (1.9 keV), 
EM = 7.1 x 10 48 cm -3 , and yf ed = 1.21 (y = 160). Right: 2-T (two-temperature) function 
(green and yellow curves with the red curve representing the sum of the two) fitted over 
the 1.5 — 8.5 keV energy range with T\ = 22 MK (1.9 keV), T 2 = 2.9 MK (0.25 keV), 
EMi = 7.6 x 10 48 cm -3 , EM 2 = 4.3 x 10 49 cm -3 , and y 9 ed = 1.18 ( v = 175). The same 
pre-flare background spectrum is shown in pink in both plots. 
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achieved over the range 1.5 — 8.5 keV, i.e., including lower energies. The values of x? ed for 
each fit are both less than 1.5, indicating a satisfactory fit to the observed spectrum in each 
case. For the 1-T fit, yf ed = 1.21 with a probability P = 3.6% of exceeding this value (165 
energy bins each 0.0375 keV wide and five free parameters giving v = 160). For the 2-T fit, 
Xred = 1-18 giving P = 5.2% (187 energy bins and 12 free parameters giving v = 175). 

The abundances of the elements were set as free parameters and determined for each of 
the functional forms given above. The Fe and Ca abundance estimates are expected to have 
higher precision (of order ±20%) than for S, Si, and Ar since Fe and Ca spectral lines form 
recognizable line features in the measured spectrum (Figured]). Separate line features from 
S, Si, and Ar are not evident with the SAX energy resolution but their contributions to the 
spectrum are considerable, as can be seen from Figure |3] Here, the contributions to the 
SAX count flux spectrum of the thermal continuum emission and the line fluxes from the 
different elements are shown with the SAX energy resolution and sensitivity. The photon 
spectrum used to calculate these different components using chianti was the 1-T model 
with a temperature of 22 MK (1.9 keV) that gave the best fit to the measured count flux 
spectrum in the 5-minute interval used for Figured] Note that the contributions from the 
different elements peak at the following energies: Fe - 1.5, 6.7, and 8 keV; Ca - 3.9 keV; S - 

2.5 keV; Si - 2.0 keV, and Ar - 3.2 keV. Other elements (notably Mg) give rise to emission 
peaking at ~ 1.6 keV. The SAX sensitivity is poorly known at such low energies because 
of the uncertain absorption of overlying material in the instrument and the presence of a 
steep and time-varying background spectrum, so spectral fits were limited to energies above 

1.5 keV and no attempt was made to determine the Mg abundance. 


As mentioned above, the continuum is made up of free-free and free-bound radiation, 
and so varying element abundances in the spectral fits in principle affects free-bound 


emission. 


Phillips & Dennis (2012) found that varying the Fe abundance from the (e.g. 
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photospheric (Asplund et ah 


20091) to coronal (IFeldman et al.lIl992l) values leads to only a 


few percent difference in the total continuum at 10 keV for a temperature of 25 MK; for the 
lower temperatures and energies mostly appropriate here, the effect will be less and can be 
ignored. 


Estimates of abundances and FIP biases of Fe, Ca, S, Si, and Ar from the two model 


fits to the spectrum of Figure [2] are give n in Tab. 


biases from the abundance sets of 


Feldman et al 


el Also given are the corresponding FIP 


(1992) and 


Fludra fe Schmelz (1999J). 


Figure [4] shows the temperatures, emission measures, and estimated abundances of 
Fe, Ca, S, Si, and Ar for all time intervals using the 2-T fit procedure over the duration 
of the 2007 June 1 flare. The temperatures and emission measures derived from the flux 
ratio of the two GOES channels are also shown. From 14:30 to 15:30 UT, there is a hotter 
component with a temperature generally above 11.6 MK (1 keV) and a cooler component 
with a temperature below 6 MK (0.5 keV). Apart from the first few time intervals of the 
flare during its impulsive phase, the GOES temperature is only slightly less than the hotter 
component (labeled “SAX Tl” in Figure EO). 


The abundance values for all elements in Figure [4] show variations apparently larger 


than the statistical uncertainties given by the OSPEX routine. T 
obtained from the error matrix and correspond to Xmd + 1 values i 


lese uncertainties are 


Bevington & Robinson 


2003). However, they do not take into account the interdepenc 


parameters, which have been recently discussed by 


Ireland et al. 


encies of the different 


(2013)), and hence should 


be considered as lower limits. The horizontal line in the plot for each element shows the 
weighted average abundance for the time interval represented by the line length when the 
count flux spectrum could be measured over the full energy range. 
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Fig. 3.— Contributions to the SAX count-rate spectrum (black curve) for the interval shown 
in Figure [2] made by the thermal continuum (red) and line complexes of elements indicated in 
the legend calculated using chianti (v. 7.0) for the 1-T model with T — 22 MK (1.91 keV) 
and with the SAX spectral resolution (598 eV at 5.9 keV). 
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Fig. 4.— Estimated emission measures, temperatures, and abundances (with ±la uncer¬ 
tainties) from SAX spectra for the 2007 June 1 flare (maximum 14:59 UT) from the two- 
temperature (2-T) fit function. Top two panels: emission measures and temperatures (keV) 
determined from SAX (red and black) and from GOES (green). Panels 3-7 from top: FIP 
biases for Fe, Ca, S, Si, and Ar. The horizonal lines show the weighted average abundance es¬ 
timates over the indicated time interval, and these estimates with uncertainties are indicated 
in the legend. 
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2.3. Comprehensive analysis of SAX Flares 2007 — 2013 


The detailed analysis of the 2007 June 1 flare above was generalized to spectra taken 
in flares detected by SAX in the 2007-2013 period and selected according to criteria given 
in Section 12.11 Care was taken to avoid selection biases and effects of varying instrument 
performance and the changing environment in which MESSENGER was located, both in 
its cruise phase and in orbit about Mercury. A semi-automatic routine was written in the 
Interactive Data Language (IDL) to analyze all the larger flares detected with SAX up 
to 2013 September 13 with the same basic techniques described using the 1-T and 2-T 
temperature functions. (The multitherm_pow and multitherm_exp functions were not used 
at this stage since acceptable fits were generally achieved with the 1-T or 2-T models.) The 
procedure handled the variations in the performance and background spectrum over this 
extensive period, analyzing spectra in each time period when the count rate was above a 
given threshold. 


The data were analyzed in three stages. First, a variable background spectrum 
was determined for each 24 hr period of the observations. This was achieved by 
selecting the 10% of all the observational time periods with the lowest count rates in 
the 1.3 — 5.3 keV energy range an d u sing those values to determine a smoothed time 


variation by applying a 


Savitzkv fe Golavl (11964]) filter with a width of 128 intervals. 


This same time profile, normalized by the average rate in the background intervals, was 
used to estimate the background flux during a flare interval for each energy bin. (See 


http://hesperia.gsfc.nasa.gov/ssw/packages/spex/doc/ospex_explanation.htm#Baekground 
for details.) The background spectrum estimated in this way was then subtracted from 
the measured spectrum in each flare time interval for analysis in OSPEX. Only those 
background-subtracted spectra having count rates in the 6.3 — 7.0 keV range of > 5000 s^ 1 
were selected for analysis, this being the limit at which the Fe-line complex at 6.7 keV could 
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be reliably measured. 

In the second stage, OSPEX was used to analyze the background-subtracted count flux 
spectrum for each selected time interval using both the 1-T and 2-T temperature models. A 
FITS file was generated for each model covering an entire one-year period of observations, 
containing the count flux spectra and best-fit model spectra with all parameters for each 
analyzed time interval. 

For the 2-T model, the spectral fit was performed over the energy range 1.5 — 8.5 keV 
as with the analysis of the 2007 June 1 flare discussed above. This range avoids both 
low energies (< 1.5 keV) with a variable instrumental component (see Figure [2]) and 
high energies with variable background count rates from charged particles. The free 
parameters include two emission measures (EM^ and EM 2 ) and temperatures (Tj and To 
where Xj < T 2 ) as well as abundances of Fe, Ca, S, Si, and Ar (in order of decreasing 
measurement certainty). For 35% of the time intervals, the emission measure of the cooler 
(Ti) component fell to values of < 10 43 cm -3 , whereupon the fitting routine automatically 
set its value to zero and continued fitting with a single temperature and emission measure. 
A pseudo-function (drm_mod) was also included to permit the automatic adjustment of the 
SAX nominal energy resolution and energy calibration. The instrument energy calibration 
can be found from the recognizable peaks at the known energies of the Fe and Ca line 
complexes, allowing both the intercept and slope of the energy calibration function to 
be estimated. The observed width of the Fe-line complex gives the detector resolution 
(FWHM, as a fraction of the measured pre-launch value). The number of free parameters 
in this case is therefore 12 (10 when EMi was set to zero). Some 186 energy bins (each 
0.0375 keV wide) are included in the 1.5 — 8.5 keV fit energy range, giving v = 174. 

For the 1-T case, we performed spectral fits over the energy range 5.5 — 8.5 keV (80 
energy bins), which is similar to the range used for the fits to RHESSI spectra done in the 
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analysis by 


Phillips & Dennis (2012). In this range, only the Fe-line complexes at 6.7 keV 


and 8 keV are included, so only the Fe abundance can be estimated (ionized Ni lines make 
a small contribution to the 8 keV complex). Also, since the Ca line complex is not included, 
only the resolution and intercept of the instrument’s energy calibration function can be 
determined. These, with T, EM, and the Fe abundance, are the free parameters, so v = 75. 


Once the FITS hies were made for all spectra and hares, selection and statistical 
analysis of the spectra were made in a third and final stage. Several criteria were used 
for including those spectra that were free from unwanted environmental, instrumental, or 
other effects. Thus, spectra were only included if: (a) they were in hare decay periods 
(those in the initial stages of hares had evidence of time-varying element abundances as 
suggested by the plots in Figure 0J; (6) for the 1-T case, temperatures were in the range 
12.8 — 23.2 MK (1.1 — 2.0 keV) (at temperatures < 12.8 MK the emission function of the Fe 
6.7 keV line complex falls sharply so the measured intensity of the line above the continuum 
is uncertain, while spectra having measured temperatures > 23 MK were nearly always the 
result of charged particle contamination of the X-ray spectrum); (c) emission measure EM 
> 10 48 cm -3 (lower emission measures gave rise to count rates too low to obtain reliable 
abundance estimates); (d) there were at least two contiguous 5-minute time intervals; (e) 
the values of Xred °f the fitted spectra were in the range 0.5 — 1.5 (values outside this 
range indicated poor fits or low count rates); (/) there was no evidence of charged particle 
contamination (identified by the count rate ratio in the 9 — 9.5 keV to 8 — 8.5 keV ranges 
being larger than 20). 


3. RESULTS OF COMPREHENSIVE ANALYSIS 

We summarize the results of the comprehensive analysis of 526 flares observed with 
SAX between 2007 May and 2013 August. The Fe abundances are presented first as 
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determined from the 1-T spectral fits over a narrow energy range, and these are followed 
by the abundances of Fe, Ca, Si, S, and Ar determined from the 2-T model fits over a 
broader energy range. (As indicated in Section 12.11 the results and plots are given in 
http://hesperia.gsfc.nasa.gov/messenger/resnlts/plots_and_tables/.) 


3.1. Fe abundance estimates using the 1-T function 

The Fe FIP bias estimates for all spectra analyzed with the 1-T function over the 
energy range 5.5 — 8.5 keV are plotted against time in Figure [5] (top panel). The period 
analyzed covered the declining portion of Cycle 23 (2007-2008) and the rising portion 
of Cycle 24 (2010-2013). With fewer free parameters, the 1-T fits are likely to be more 
reliable than those from the 2-T function (Section 13.2ft . Plotted in Figure [6] (top panel) 
are FIP bias values plotted against flare sequence number, with the FIP bias averaged 
over the duration of each flare. This plot shows the existence of flare-to-flare variations 
since the scatter of points is significantly larger than the uncertainties on individual points. 
The uncertainties are a combination of the uncertainties given by the OSPEX routine 
and variations between different times during the same flare. There is at least one group 
of flares, flare numbers 238-261 occurring between 2011 October 28 and 30, that have 
mostly smaller-than-average Fe abundances. These flares may all come from a single active 
region although MESSENGER was then viewing the solar hemisphere partly directed away 
from the Earth so we cannot be sure of this. STEREO B was viewing the Sun from a 
similar angle to MESSENGER at that time and showed an active region at a Stonyhurst 
hcliographic longitude of ~ —100° that could have been the origin of at least some of these 
flares. None of them are evident in the GOES light curves except for six on 2011 October 
28 all at the Cl to C2 level. 

Figure [7] (top panel) shows the frequency distribution of FIP bias estimates obtained 
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Fig. 5.— Top: Fe FIP bias plotted against time estimated from spectral fits with the 1-T 
function in the energy range 5.5 — 8.5 keV. Bottom: FIP bias over the same period from 
spectral fits with the 2-T function in the energy range 1.5 — 8.5 keV. The broken red line at 
the top shows the SAX data coverage. The green and blue ho ri zonta l lines show the mea n 


FIP bias with ±lcr uncertainties given by Phillips & Dennis (2012) and Warre n (2014). 
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Fig. 6.— Top: Fe FIP bias determined from 1-T function fits (fit range 5.5 — 8.5 keV) 
averaged over each flare plotted against flare sequence number. Th e gree n and blue horiz ontal 


lines show the mean FIP bias and ±lcr uncertainties given by 


Phillips & Denn is (2012) and 


Warren] (20141 respectively. Middle: same as the top plot with 2-T function fits (fit range 


1.5 — 8.5 keV) with the abundances of other elements set free. Bottom: same as the middle 
plot but for the Ca FIP bias. 
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for all spectra. From this distribution, the mean FIP bias from the 1-T function is 
1.62 ± 0.28 (s.d.), corresponding to an Fe abundance (on a logarithmic scale with H = 12) 
of A(Fe) = 7.71 ± 0.07. The standard deviations include any real time variations over the 
period of observations. 

To check for any trends with flare parameters, the Fe FIP bias is plotted against both 
emission measure and temperature in Figure HI There appears to be a slight tendency 
for lower FIP biases to result from larger emission measures (> 10 49 cm -3 ) and lower 
temperatures, but the significance of the trend is low. 

The occurrence distribution of Xi- e d va l ues is of interest as it checks whether there is 
a match with that expected for a chi-squared distribution with the appropriate number 
of degrees of freedom {y = 75: Section ITTSjh Figure 0 (top panel) shows the observed 
distribution for the 1-T fitting function and the statistically expected distribution for 
v = 75. There is a very close agreement of the observed and expected distributions, showing 
that lcr-normalized residuals for each energy bin are close to expected and that there are no 
significant instrumental effects in the 5.5 — 8.5 keV energy range affecting the distributions. 


3.2. Abundance estimates using the 2-T function 

Abundance estimates of Fe, Ca, S, Si, and Ar were made using the 2-T fitting function 
over the broader energy range of 1.5 — 8.5 keV. The three free parameters describing the 
SAX detector energy resolution and calibration were found to vary by less than 20% over 
the entire 7 yr period of the observations, indicating the stability of the SAX instrument. 

Estimates of the Fe FIP bias from the 2-T analysis are summarized in Figure 0 
(bottom panel) and Figure [6] (middle panel). Figure [7] (bottom panel) shows the frequency 
histograms for Fe FIP bias values from this analysis. While there is little difference in 
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Fig. 7.— Frequency histograms for the Fe FIP bias values determined from fits using the 
1-T function (top: energy range 5.5 — 8.5 keV) and 2-T function (bottom: 1.5 — 8.5 keV) for 
all selected time intervals. The mean FIP bias is 1.62 ± 0.28 (s.d.) for the 1-T distribution 
and 1.66 ± 0.50 for the 2-T distribution. 
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Fig. 8.— Fe FIP bias estimates from the 1-T fitting function (energy range 5.5—8.5 keV) from 
all spectra plotted against emission measure (top) and against temperature (bottom). In 


both plots, th e g reen a nd b lue horizonta 


lines show the mean FIP bias with ±1 a uncertainties 


given hv IPhillips fc Dennisl (120121 1 and IWarrenl (120141 1. respectively. 
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Fig. 9.— Top: Frequency distribution (shown in black) of the reduced chi-squared values 
for all time intervals fitted with the 1-T function (energy range 5.5 — 8.5 keV) with expected 
distribution (in red) for v = 75. Bottom: Frequency distribution of reduced chi-squared 
values for all time intervals using the 2-T function (energy range 1.5 —8.5 keV), with expected 
distribution (in red) for v = 174. 
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the mean value of FIP bias, 1.66 ± 0.50 (s.d.), or A(Fe) = 7.72 ± 0.11, the amount of 
scatter is larger. As with the 1-T case, there is little apparent dependence of Fe FIP bias 
on emission measure or temperature of the higher-temperature component. The frequency 
distribution of the reduced chi-squared values obtained for all analyzed time intervals in 
the 2-T case deviates significantly from the expected distribution for v = 174 (Figure [9]). 
This difference did not show up for the 1-T fits probably because of the narrower energy 
range. It most likely reflects some small instrumental effect not included in the instrument 
response matrix. The values of the lcr-normalized residuals averaged over all analyzed 
spectra for each energy bin show that certain bins consistently have larger or smaller count 
fluxes at the few percent level but the origin of these small discrepancies is unknown. Since 
the difference is only ~ 0.1 in the reduced chi-squared value, it is does not result in any 
significant deviation in the mean values of the free parameters obtained from the fits. 

Abundance and frequency distribution plots were also obtained from the 2-T function 
for Ca, S. Si, and Ar. Figure [6] (bottom panel) shows the mean Ca FIP bias for each 
flare plotted against flare sequence number. There is a slight but statistically significant 
decrease for events 238-261, as with the Fe FIP bias values shown in this plot (top two 
panels). Frequency histograms for Ca and also for S, Si, and Ar are given in Figure [10J The 
mean Ca FIP bias is 3.89 ± 0.76 corresponding to A(Ca) = 6.9 ±0.1 with the uncertainty 
including any possible real variations from flare to flare. 

Figure [10] (top right and bottom panels) show frequency histograms for S, Si, and Ar. 
As these elements do not have recognizable line complexes in SAX spectra, the FIP bias 
determinations are less reliable. We obtained mean values for all flares to be 1.23 ±0.45 (S), 
1.64 ± 0.66 (Si), and 2.48 ± 0.90 (Ar). These correspond to abundances (on a logarithmic 
scale with H = 12) of A(S) = 7.2 ± 0.2, A(Si) = 7.7 ± 0.2, A(Ar) = 6.8 ± 0.2. 

The mean FIP bias values and corresponding element abundances determined from 


Number of Occurrences per bin 


Ca FIP Bias 


S FIP Bias 



Si FIP Bias Ar FIP Bias 


Fig. 10.— Frequency histograms for the FIP bias values for Ca (top left), S (top right), Si 
(bottom left), and Ar (bottom right) determined from the 2-T fitting function (energy range 
1.5 -8.5 keV). 
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this analysis are given in Table 2. The FIP bias is plotted against the FIP for each element 
in Figure |TT] 


DISCUSSION AND CONCLUSIONS 


The FIP bias results from our analysis of MESSENGER SAX solar flare data cover the 
period 2007-2013, when the spacecraft went from the vicinity of the Earth to orbit about 
Mercury. They are summarized in Table [21 

Other investigators have estimated flare abundances of the elements considered here 
by various methods and so it is of interest to compare their results with ours. RHESSI 
is a broad-band X-ray spectrometer like SAX covering the energy range above ~5 keV 
with a slightly poorer spectral resolution of ~ 1 k eV but still adequate to resolve the 
Fe-linc complex at 6.7 keV. Analysis by Phillips & Dennis ( 2012 ) of the decay phase of 


20 strong flares (GOES class up to X8) seen with RHESSI gave A(Fe) = 7.91 ± 0.10 or 
a FIP bias of 2.6 ± 0.6. This is higher than the result obtained here though just within 
the statistical uncertainties of each determination. The X-ray Solar Monitor (XSM) on 
the Indian lunar space mission Chandrayaan-1 observed solar 1.8 — 20 keV X-ray emission 
to complement the X-ray fluorescence instrument on board designed to investigate the 
lunar surface composition. XSM had a spectral resolution of 200 eV, an improvement 
over SAX and enough to resolve spectral features of all the elements reported here. The 
mission only worked for nine months during a solar-quiet period and n o str ong f lares were 


observed. However, a C2.8 flare in 2009 was analyzed by 


Narendranath et ah 


(2014), 


and they reported a weighted avera ge A (Fe) = 7.99 ± 0.09 corresponding to a FIP bias 


of 3.1 ± 0.7, which is similar to the 


Phillips & Dennis (2012) RHESSI resu lt but higher 


than the SAX value obtained here. For several smaller flares, 


Narendranath et al. 


(2014) 


reported A(Ca) = 6.6 ± 0.2 and A(S) = 7.3 ± 0.3, consistent with our abundance estimates. 
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Table 1. Abundances and FIP biases from the SAX spectrum of the 2007 June 1 flare 

(14:56:30-15:01:29 UT) 


Element FIP Fitting Model (see text) Photo- Corona 6 



(eV) 

1-T 

2-T 

aT~ a 

a exp(— T/t) 

sphere 0. 


Fe 

7.87 

2.33 c [7.87 d ] 

1.70 [7.73] 

1.22 [7.59] 

1.76 [7.74] 

7.50±0.04 

3.98 [8.10] 

Ca 

6.11 

3.89 [6.93] 

3.59 [6.90] 

3.17 [6.84] 

2.98 [6.81] 

6.34±0.04 

3.89 [6.93] 

S 

10.36 

0.96 [7.10] 

0.78 [7.01] 

0.46 [6.78] 

0.45 [6.77] 

7.12±0.03 

1.41 [7.27] 

Si 

8.15 

2.04 [7.82] 

1.09 [7.55] 

0.69 [7.35] 

0.75 [7.39] 

7.51±0.03 

3.89 [8.10] 

Ar 

15.76 

0.85 [6.33] 

1.48 [6.57] 

0.71 [6.25] 

0.37 [5.97] 

6.40±0.13 

1.51 [6.58] 


a Photospheric abundances from 


Aspl un d et a l. 


(2009) on logarithmic scale with A(H) = 12. 


b CHIANTI “coronal abundances” from 


Feldman et al 


( 19921) . 


C FIP bias with uncertainties (abundances on logarithmic scale in square brackets). 
d Abundances on a logarithmic scale with A(H) = 12. 


Table 2. FIP Biases from SAX spectra for 526 flares (2007 May 28-2013 August 19). 



FIP 

(eV) 

Photo- 0 

sphere 

CHIANTI 

Corona Hybrid 

RESIK 

RHESSI 

SAX 

2007 d 

SAX 

2007-13 

XSM 6 

SDO/EVE c 

Fe 

7.87 

7.50 

3.98 

2.14 

2.57 

1.70 

1.66±0.34 

3.09=t0.80 

1.27±0.22 

Ca 

6.11 

6.34 

3.89 

2.09 

3.89 

3.59 

3.89±0.76 

1.82±1.06 


S 

10.36 

7.12 

1.41 

1.58 

1.10 

0.78 

1.23±0.45 

1.51±1.51 


Si 

8.15 

7.51 

3.89 

2.29 

2.34 

1.09 

1.64±0.66 

2.63±1.64 


Ar 

15.76 

6.40 

1.51 

0.91 

1.10 

1.48 

2.48±0.90 



K 

4.34 

5.03 

4.37 

2.69 

6.76 






a Absolute abundances from 


Asplund et al, 


( 20091 ) on logarithmic scale with A(H) = 12. 


b FIP bias values from 
C FIP bias values from 


Narendranath et al 


( 20141) . 


Warrenl ( 20141) . 


d FIP bias values from Table [I] for the 2-T fits to the 2007 June 1 flare. 
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Fig. 11.— Measured FIP bias values as a function of FIP. Red circles: current results for Fe, 
Ca, Si, S, and Ar (FIP bias values moved to the right by 0.15 eV for clarity). The error bars 
show the ±lcr values obtained from the scatter of the measured values shown in Figure [7] for 


Fe and Figure [TO] for other elements. Also shown are the following pu 


black plus signs: CHIANTI “coronal” abundance set 


CHIANTI “hybrid” a 
SDO/EVE given by 


umdance set 


Fludra fe Schmelz 


(Asplund et al. 


dished FIP bias values: 


2009[); green asterisks: 


19991 1: purple square: Fe FIP bias from 
Warren l (20141 (moved to the left by 0.15 eV for clarity); blue p lus signs 


with ±lcr error bars: RES IK K, Ca, Si, S, Ar valu es given by 


RHESSI Fe value given by 


S vlwester et ah 


(2014) and 


Phillips & Dennis (2012). 
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Broad-band X-ray spectrometers, in general, have an inferior spectral resolution to 
that of crystal spectrometers flown on numerous spacecraft since the 1960s. However, 
they have the advantage that the continuum is unambiguously observed whereas many 
crystal spectrometers suffered from a large background formed by the fluorescence of 
the crystal material by solar X-rays. This difficulty has been overcome in some cases 
by using crystals made with material having low atomic number, as with the graphite 
crystals in the OSO-8 spect rometer. Solar flare data from this instrument were analyzed 


by 


Veck & Par kin son ( 1981 1 giving A(Si) = 7.7!°!, A(S) = 6.9ztoil, A(Ar) = 6.4to;3, and 


A(Ca) = 6.5io;2- While the Si abundance is very close to our value, S, Ar, and Ca are all 
lower than ours and just outside the range indicated by the error bars. The fluorescence 
background was reduced to zero or to an accurately measured amount in the RESIK crystal 
spectrometer. A differential emission measure analysis of RESIK observations of a large 


flare gave A (Si 
which, like the 


= 7.56 ± 0.08, w 


rich is consistent with our result, and A(S) = 6.94 ± 0.06 


Veck & Parkinson! (119811 ) result, is lower than our result. Earlier analysis 


of the Ca abundance during many flares observed by the Solar Maximum Mission Bent 


Crystal Spectrometer by 


Sylwester et al. 


(1998) showed that there was a factor-of-three 


variation from flare to flare, with a mean value A(Ca) = 6.76 ± 0.10, consistent with our Ca 
abundance and its uncertainty. Thus, in the light of this result, the apparent flare-to-flare 
abundance variations found here (e.g., the slight decrease for flares 238-261) may be real. 


Warren (120141 ) used the fluxes of lines due to Fe ionization stages from Fe XV to 


Fe XXIV obtained from the Solar Dynamics Observatory Extreme ultraviolet Variability 
Experiment (EVE) (wavelength range: 60-200 A). His result for the FIP bias of Fe using 


a DE 
or by 


VI analysis was 1. 


7± 0.22, corresponding to A(Fe) = 7.57, lower than obtained here 


Phillips & Dennis (j2012h . The form for the temperature dependence of the DEM was 


taken to be the sum of Gaussian functions fixed at particular values of log T, an artificial 
assumption which is far from the near-bimodal functions derived by other authors for 
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solar X-ray flares (e.g., iMcTiernan et al.lll999t ISvlwester et al.ll2014j) . Also, there is little 
constraint on the high-temperature part of the DEM function which may lead to some 
uncertainty in its determination, while the crowded nature of emission lines in the (e.g.,) 
100 — 120 A region is only just resolved by the ~ 1 A spectral resolution of EVE. 

Elemen t enhanc e ment s in closed coronal loops by a FIP effect have been theoretically 


modeled by 


Laming (120041) . The idea is that the more easily ionized low-FIP elements 


are preferentially accelerated into the corona by a ponderomotive force associated with 
Alfven waves that travel in the loop down to the chromospheric footpoints. The latest 


evolution of this theory also includes slow-mode waves ( Lamin 


2012). The particular case 


in these calculations most appropriate to solar flares indicates enhancements by factors 
of 2-3 (relative to oxygen) for Ca and Fe and more modest enhancements for S and Si 
and no enhancement for the high-FIP element Ar. Apart from Fe for which the predicted 
value is larger than obtained here, our results have some resemblance to the predictions of 
one of the model calculations with small wave amplitudes. However, the 100,000 km loop 
length in the calculations is much larger than typical X-ray flare dimensions as revealed 
by RHESSI images, so that a more detailed comparison with realistic loop dimensions and 
other parameters is clearly desirable. 

The small enhancements over photospheric abundance indicated in this work for the 


relatively high-FIP e leme nts S and Ar are con s isten 


Feldman et al. 

(1992 

) and 

Feldman & Laming 

(2000) 


; with the original ideas expressed by 


but our low enhancements for Si and Fe (FIP values ~ 8 eV) are not expected. The much 
higher enhancement that we find for Ca (FIP = 6.1 eV) still indicates the presence of a FIP 
effect for solar flare plasm a, as does the very high enhancement factor found from RESIK 


spectra of 5.8 ±1.6 ( Svlwester et al 


20101) ) for K with its very low FIP of 4.3 eV. This would 


seem to point to a boundary between low-FIP and high-FIP elements that is significantly 
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less than the 10 eV that is normally assumed, perhaps about 7 eV. This would imply that 
the fractionation takes place at a l ower te m peratu re than has been conceived of in the past. 


Results for a white-light flare ([Martinez Oliveros et al 


2012) point to a similar low altitude 


in the solar atmosphere of white-light emission and a hard X-ray (30-40 keV) source as 
observed by RHESSI for a near-limb flare, around 200-300 km above the photosphere, 
below the temperature minimum and considerably below the altitude expected from the 
standard thick-target model for hard X-ray emission. If this applies to flares generally, the 
origin of the flare plasma may be at lower chromospheric altitudes than had been previously 


thought. Likewise, the deviations from photospheric abunda nces in fl are p 


from a fractionation process, perhaps like that described by 


asrnas may arise 


Lamind (120041) but operating 


at lower temperatures than has previously been considered. The apparently time-varying 
Ca and Fe abundances from flare to flare (as indicated for both Ca and Fe in this work and 
Ca from 


Sylwester et ah 


(119981 )) should also be considered in any theoretical modeling of 


the FIP effect. 


We acknowledge the MESSENGER XRS project team for making all data promptly 
and publicly available for independent analysis. CHIANTI is a collaborative project involving 
the US Naval Research Laboratory, the Universities of Florence (Italy) and Cambridge 
(UK), and George Mason University (USA). 
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